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Abstract 

We report the result of a systematic methanol observation toward IRAS 
19312+1950. The properties of the SiO, H 2 O and OH masers of this object are con¬ 
sistent with those of mass-losing evolved stars, but some other properties are difficult 
to explain in the standard scheme of stellar evolution in its late stage. Interestingly, a 
tentative detection of radio methanol lines was suggested toward this object by a pre¬ 
vious observation. To date, there are no conhrmed detections of methanol emission 
towards evolved stars, so investigation of this possible detection is important to better 
understand the circumstellar physical/chemical environment of IRAS 19312+1950. In 
this study, we systematically observed multiple methanol lines of IRAS 19312+1950 
in the A = 3 mm, 7 mm, and 13 mm bands, and detected 6 lines including 4 ther¬ 
mal lines and 2 class I maser lines. We derived basic physical parameters including 
kinetic temperature and relative abundances by htting a radiative transfer model. 
According to the derived excitation temperature and line prohles, a spherically ex¬ 
panding outflow lying at the center of the nebulosity is excluded from the possibilities 
for methanol emission regions. The detection of class I methanol maser emission sug¬ 
gests that a shock region is involved in the system of IRAS 19312+1950. If the central 
star of IRAS 19312+1950 is an evolved star as suggested in the past, the class I maser 
detected in the present observation is the hrst case detected in an interaction region 
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between an evolved star outflow and ambient molecular gas. 

Key words: maser — stars: individual (IRAS 19312+1950) — stars: mass loss 
— radio lines: stars 


1. Introduction 

SiO maser sources are, in most cases, identifled as being associated with mass-losing 
evolved stars of either asymptotic giant branch (AGB) stars, post-AGB stars or red super¬ 
giants (RSGs). In fact, the maser properties of SiO, H 2 O and OH masers associated with 
IRAS 19312+1950 (119312, hereafter) are consistent with those of AGB/post-AGB stars or 
RSGs (Nakashima et ah 2011), and the properties of the central infrared point source are also 
consistent with those of mass-losing evolved stars (Murakawa et ah 2007). Nevertheless, some 
other characteristics of this object cannot be explained in the standard scheme of late-stage 
stellar evolution. For example, an extended nebulosity with point symmetric structure, which 
shows complicated knots, surrounds a spherically expanding molecular outflow lying at the cen¬ 
ter (Nakashima & Deguchi 2000; Nakashima & Deguchi 2004; Deguchi et ah 2004; Nakashima 
& Deguchi 2005; Murakawa et al. 2007), the gas-mass derived from CO lines is relatively high 
for an evolved star (10-15 Mq; Deguchi et al. 2004), and a rich set of molecular species has 
been detected toward the nebulosity (so far 22 molecular species have been detected; Deguchi 
et al. 2004; Nakashima & Deguchi 2005). 

Since SiO maser emission is occasionally detected toward dense cores in molecular clouds 
(for example, Ori IRc 2, W51 IRs 2, and Sgr B2 MD5; Hasegawa et al. 1986), which are 
detectable in various molecular lines, one possible explanation might be that 119312 is a young 
stellar object (YSO) embedded in a molecular cloud. However, the characteristics of 119312 
are clearly different from such SiO maser sources embedded in molecular clouds. For instance, 
119312 shows an isolated point-like feature in mid-infrared images (see, e.g., Nakashima et al. 
2011), while all the SiO maser sources found in molecular clouds are embedded in a huge star 
forming region, of which the background is generally very bright at mid-infrared wavelengths. 
Additionally, there are no clear star forming activities in the vicinity of 119312, such as an 
enhancement of the number density of stars. 

As of now, three hypotheses are being suggested as possible explanations for the origin 
of 119312 (see, e.g., Nakashima et al. 2011): (1) a red nova formed by the merger of two main 
sequence stars (or two stars going to main sequence), (2) a mass-losing evolved star physically 
associated with a small, isolated molecular cloud, and (3) a mass-losing evolved star behind a 
small isolated molecular cloud, which is not physically associated with the star (but the star 
and molecular cloud shares the same systemic velocity). However, none of these hypotheses 
can adequately explain the observational characteristics of 119312. 


2 


Interestingly, Deguchi et al. (2004) reported a tentative detection of radio methanol lines, 
bnt this tentative detection is unexpected because the most likely interpretation of 119312 are 
that it is a mass-losing evolved star. In fact, there have been no previous secure detections of 
radio methanol lines toward any evolved stars, though some sensitive methanol searches have 
been made toward the samples of evolved stars (see, e.g., Charnley & Latter 1997; Ford et ah 
2004; He et al. 2008). 

If the methanol detection were to be conhrmed it would be of great signihcance in studies 
of the physical and chemical environments of evolved stars. The detection of multiple methanol 
lines enables essential physical parameters of the emission region to be determined, which would 
be invaluable in understanding this unusual source. 

In the present research, we systematically observed multiple methanol lines at A = 3 mm, 
7 mm, and 13 mm bands toward 119312, and securely detected 6 methanol lines including 4 
thermal lines and 2 class I maser lines. The detected class I masers seem not to be classihed 
into any known categories of methanol masers. We obtained some basic physical parameters 
by constructing a simple radiative transfer model. Based on the results, we discuss the nature 
of the methanol emission region, implying possible locations of the emission region of methanol 
lines. 

2. Details of Observation and Results 

The observations were made with the Nobeyama Radio Observatory (NRO) 45 m tele¬ 
scope on 2012 March 14-17 and May 14-16. [Only the 7(0,7)—6(1,6) A++ line was separately 
detected on 2011 February 11.] The coordinate values of the target, which were used in the 
observation, were 19^33™24®249, -|-19°56'55”65 (J2000.0). The receivers used in the present 
observations are summarized in table 1. Since the frequency coverage of H40 (2 GHz) is much 
wider than S40 (0.5 GHz), in some cases we used H40 to cover multiple lines at a time to shorten 
the total observing time, even though the system temperature of H40 is slightly higher than 
that of S40. The half-power beam widths (HPBW) of the telescope at A =13 mm, 7 mm and 
3 mm were roughly 80”0 (at 22 GHz), 39"6 (at 43 GHz) and 16"l (at 110 GHz), respectively. 

The acousto optical spectrometers, AOS-H and AOS-W, were used, and the details are 
summarized in table 2. Since the line intensity of the observed methanol lines was relatively 
weak, we mainly used AOS-W, of which the velocity resolution is lower than AOS-H. The 
antenna temperature (T^) used in the present paper is corrected for the atmospheric and tele¬ 
scope ohmic loss, but not for the beam nor aperture efficiency. Observations were made in a 
position-switching mode, and the off-position was chosen to be a location -|-20' away from the 
object in the right ascension, where contamination is minimized (the background emission at 
the off-position was carefully examined in our previous study; see, Nakashima et al. 2004). Data 
reduction, including flagging bad channels and subtraction of the baseline, was made using the 
software package NEWSTAR, which has been developed by NRO. 
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In figure 1, we present the spectra of 6 detected lines, and the line parameters of the de¬ 
tected lines are summarized in table 3; the line parameters include rest frequencies, transitions, 
intensity peak velocities, peak intensity, velocity integrated intensities. The rms noise (in T*) 
and on-source integration time are also given in table 3. For non-detections, rms noises (in T*) 
are are given in table 4. We use the notation of transitions used in Lovas (1992) throughout 
the paper. The detected 6 lines include 4 thermal lines (48.37 GHz, 96.739 GHz, 96.741 GHz, 
and 96.745 GHz) and 2 maser lines (44.07 GHz and 95.17 GHz). The two detected maser 
lines belong to class I methanol masers (see, e.g., Menten 1991). As we mentioned above, 
the 44.07 GHz maser line was separately observed on 2011 February 26, which was roughly 9 
months before the observation of the other lines. The integration time and achieved rms noises 
of the observation on 2011 February 26 is given in table 3, and those on 2012 May are 940 sec 
and 0.013 K (i.e., the later observation on 2012 May was more sensitive than the previous one 
on 2011 February). We observed this maser line again in 2012 May, but the result was negative. 
In addition, we observed 2 class H methanol lines (i.e., 23.121024 GHz and 108.893929 GHz; 
Wilson et ah 1984; Voronkov et al. 2011), but again the results are negative in the present 
observation. 

All 4 thermal lines exhibit a linewidth of about 3-4 km s“^, while the maser lines 
show a narrower linewidth of about 1-2 km s“^ (here, the linewidths are measured at the 0 - 
intensity level). The intensity peak velocities of detected lines exhibit consistent values of about 
35.5 km s“^ - 36.1 km s“^ except for the 96.744 GHz line, of which the intensity peak velocity 
is 39.1 km s“^. This small shift in the velocity is due presumably to the weak intensity of this 
line (i.e., low signal-to-noise ratio). In hgure 1, the intensities of the 1(0,1)—0(0,0) A^^ and 
8(0,8)—7(1,7) A+’'' lines are relatively weak compared with other lines (signal-to-noise ratios are 
roughly 3-4). However, the intensity peak velocities of these lines almost exactly correspond to 
those of other lines, and additionally the emission is detected in multiple frequency channels; 
this situation supports that the emission is originated in a natural astrophysical object rather 
than spurious emission, which is often detected in a single frequency channel. One may think 
that these could be spurious emission, by chance, lying exactly at the systemic velocity of the 
object. However, the chance probability of hnding such very unusual artihcial features in two 
frequency bands at a time would be negligibly small. Thus, it is natural to conclude that these 
two lines are the true emission of methanol molecules. 

In the bottom panel of hgure 1, we see another possible emission at about —25 km s“^. 
This emission exhibits the intensity of 0.032 K with a signal-to-noise ratio of 4.6 a. This 
emission is close to the G 4 H 10-9 J=19/2-17/2 line at 95.18894 GHz (the radial velocity of the 
emission is 36.5 km s“^ if we assume the rest frequency of the G4H 10-9 J=19/2-17/2 line). 
However, we harbor doubts on this marginal detections for the following two reasons. Firstly, 
we did not detect another G4H line (10-9 J=21/2-19/2 at 95.15032 GHz), which is within the 
frequency range of the bottom panel of hgure 1. In fact, the G4H 10-9 J=21/2-19/2 line was 
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detected toward IRC+10216 with almost the same intensity with the C 4 H 10-9 J=19/2-17/2 
line (Gnelin et al. 1978). Secondarily, the linewidth seems to be too narrow as a thermal line; 
this nnnsnally narrow linewidth is reminiscent of the characteristic of an artihcial spnrious 
emission (this is obvious when we compare with the line prohles of other thermal lines; see, 
Appendix 1, for example). Thus, as the detection is quite marginal, we do not go deeper into 
scientihc discussions about this emission. 

3. Radiative Transfer Modeling of Detected Methanol Lines 

For better understanding of the emission region of methanol lines, we constructed a 
simple radiative transfer model, which uses the large velocity gradient (LVG) approximation. 
Dust emission and absorption within the emission region was taken into account in the way 
described in Sutton et al. (2004). We assumed that the dust particles are intermixed with gas 
molecules. The same physical temperature of the gas and dust components is assumed. The 
molecular emission region was assumed to be spherically symmetric (in the optical sense) and 
uniform in H 2 density, gas and dust temperature, gas-to-dust ratio and methanol fractional 
abundance. The influence of external infrared sources was not considered. The dust opacity 
law was chosen as Tdust ~ We adopted a gas-to-dust mass ratio of 100 and a cross section 
at 1 mm of 2.6 x 10“^® cm^ (Sherwood et al. 1980). 

In addition to the model described in Sutton et al. (2004), we used collision transition 
rates based on the model of collisions of methanol molecules with He and para-H 2 molecules 
(Gragg et al. 2005) in this work. Recently, newer collisional rates have been published (Rabli & 
Flower 2010a; Rabli & Flower 2010b), but we would like to note that the difference caused by 
using the newer rates is negligible for low transition methanol lines, such as the lines analyzed 
in the present research. The scheme of energy levels in this model includes rotational levels 
with quantum numbers J up to 22, \K\ up to 9; the levels include the rotational levels of the 
ground, hrst and second torsionally excited states. In total, 861 levels of A-methanol and 852 
levels of E-methanol were considered according to Gragg et al. (2005). 

Through the course of model htting, we made estimates of the values of hydrogen number 
density, uhj [cm“^], specihc column density of methanol, A^chsOh/^R [cm“^s], gas kinetic 
temperature, Tk [K], using the measured values of “quasi-thermal” (i.e. not maser) methanol 
lines. The brightness temperatures of three lines from quartet lines at 96.741 GHz-96.744 GHz 
and 48.37 GHz, as well as upper limits for the brightness temperatures of 10 other lines in the 
frequency range 36.169 GHz-108.893 GHz were considered (see, tables 5 and 6 ). 

In order to estimate these physical parameters, we have searched for a set of parameters 
which exhibits the best agreement between the values of the calculated brightness temperatures 
(T™°d) and the measured brightness temperatures (T°^®). This corresponds to Ending the min¬ 
imum of — T“°d)/(Tj)^, where cij is the observational uncertainty for a particular 

line. 
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The process was done in two steps. In the hrst step we explored the parameter space 
in order to hnd the approximate location of the minimnm. For this pnrpose, we used the 
database of population numbers for the quantum levels of methanol described in Salii (2006). 
We give a brief description of this database below. 

The kinetic temperature in the database ranges from 10 K to 220 K with a 10 K step. The 
hydrogen number density varies from 10^ to 10® cm“^ with a step of 0.5 dex. The fractional 
abundance of methanol molecules relative to molecular hydrogen in the population number 
database varies from 10“® to 10“® with a 1 dex step. The lower values of these parameters 
correspond to the physical conditions of dark molecular clouds, while the higher values of both 
parameters can occur in shocked molecular material. The outer parts of AGB star envelopes 
have values of kinetic temperature and hydrogen number density corresponding to the higher 
ends of the given intervals. The values of the methanol specihc column density (A^chsOh/^I^ 
[cm“^s]) in the database range from 10® cm“^ s to 10^® cm“^ s with a step of 0.25 dex. 

The variable parameter of the source size is introduced in order to take into account 
beam-dilution effects. According to the results of CO mapping observations (Nakashima et al. 
2004; Nakashima & Deguchi 2005), the source in molecular line emission seems to have a size 
of the order of 5" - 15". The analysis made after the hrst step of the search for the best ht 
parameters provided the following constraints for the model parameters: kinetic temperature 
Tk not higher than 70 K, uh is about 10^ cm“®, specihc methanol column density is about 10®, 
and the fractional abundance of methanol is not less than 10“^, the source size is about 8.6". 

At the second step of the search we made more accurate estimates of the physical pa¬ 
rameters with the constraints that were obtained in the hrst step. We found that within these 
constraints the ht is not sensitive to fractional abundance values within the interval from 10“^ 
to 10“®. For the analysis, we adopted the value for fractional abundance of 10“^. We also 
hxed the value for the source size at 8.6". Since the expected values of the kinetic temperature 
are rather low (T^ not higher than 70 K), we reduced the scheme of the energy levels in or¬ 
der to increase the robustness of the model ht and decrease the computing time. The scheme 
of methanol energy levels used in the second step was restricted to rotational levels with the 
quantum numbers J up to 18, \K\ up to 8; the scheme included rotational levels of the ground 
and the hrst torsionally excited states. In total, 380 levels of A-methanol and 374 levels of 
E-methanol were considered. 

Consequently, the set of parameters = 36 K, A^^chsOh/^I^ = 2.24 x 10® cm“^s 
(A^chsOH = 3.83 X 10^® cm“^), = 3.98 x lO'^ cm“® was found to provide the best ht with 

= 0.08. Upper limits on for undetected lines are within the error bars [it is noteworthy 
that the 23.445 CHz, 36.169 CHz, 44.069 CHz, 95.170 CHz lines experience population inver¬ 
sion (maser ehect) with this set of parameters; see, table 6]. To estimate the conhdence level of 
the ht, we used the method given by Lampton et al. (1976). According to Lampton et al. (1976) 
for 3 free parameters, a la conhdence interval includes parameters within contours + 3.5. 
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Thus, la confidence intervals are for Tk = 20-40 K, iVcHsOH/Al/ = (1.9-3.3) x 10^ cm“^s, urj = 
(0.05-2.5) X10^ cm“^. Estimated physical parameters are typical for molecular clumps experi¬ 
encing the influence of moderately strong shock waves. The results of the present calculation 
are visually summarized in figure 2. In table 6 , the theoretically predicted intensity of the 
0(0,0) —1(—1,1) E line (108.8 GHz) is 3.5 times higher than the observational 3 a rms level. 
This inconsistency is due presumably to the simpleness of the present model. Since the optical 
depth of this line is significantly larger than other lines, the intensity must be sensitive to the 
structure of the emission source. In addition, for higher frequency lines, the contribution from 
dust emission, which is not considered in the present calculation, would be non-negligible, and 
the dust emission could be stronger in higher frequencies than in lower frequencies. Therefore, 
to explain all the lines including higher frequency lines, in future we need to consider additional 
physical conditions, such as the structure of the emission sources and the emission of dust with 
the dependency on frequencies. 

4. Discussion 

In the present observations, we securely detected methanol lines toward 119312, and 
derived some informative results: e.g., kinetic temperature, line profiles and detections of class 
I maser lines, etc. In the first subsection we compare the sensitivities between the present 
observations and the previous methanol observations of evolved stars, because many char¬ 
acteristics of 119312 are consistent with those of mass-losing evolved stars. This process is 
important, because there have been no detections of methanol emission toward evolved stars. 
We need to clarify the peculiarity of 119312 by comparing the present observation with the 
previous methanol observations of evolved stars. In the next subsection, we discuss the nature 
of methanol emission and the possible locations of the emission regions based on the derived 
results. 

4.1. Comparison of Sensitivities 

As we previously mentioned in Section 1, Charnley & Latter (1997) searched the CH 3 OH 
2(0, 2)—1(0, 1) A++ line at 96.741377 GHz toward 6 evolved stars (i.e., IRG-l-10420, R Gas, 
IRG-l-10011, TX Gam, IK Tan, and OH 231.8-1-4.2). Since we also observed this line, we 
shall compare the sensitivities. The rms noise level achieved by Gharnley & Latter (1997) 
was typically 0.0028 K. If we assume a conversion factor to Jy is 35.0 S'y[Jy]/T^[K] (see, e.g., 
Papadopoulos et al. 2012) for the observation of Gharnley & Latter (1997), their typical rms 
converts to 0.098 Jy. On the other hand, the rms noise level of the present observation in 
the GH 3 OH 2(0,2) —1(0,1) A++ line is 0.009 K (T^), and this is converted to 0.036 Jy (the 
conversion factor to Jy is given in table 3). This means that the present observation is roughly 
3 times more sensitive than Gharnley & Latter (1997). However, we should note that even if 
119312 was observed in Gharnley & Latter (1997), the signal-to-noise ratio could be still more 
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than 3cr (strictly, 3.la); therefore, the emission of 119312, could, in principle have been detected 
by Charnley & Latter (1997) if they had observed it. 

Charnley & Latter (1997) searched for CH3OH 2(0, 2) —1(0, 1) A++ line emission toward 
a proto-planetary nebulae, OH 231.8-1-4.2 (Rotten Egg Nebula), which in our previous work we 
identihed as having some similarities with 119312 (Nakashima & Deguchi 2000; Nakashima & 
Deguchi 2003; Nakashima & Deguchi 2004; Nakashima et ah 2004; Nakashima & Deguchi 2005; 
Nakashima et al. 2011). The 12 micron IRAS flux density of OH 231.8-1-4.2 (18.98 Jy) is similar 
to that of 119312 (22.47 Jy), and both exhibit very red mid-infrared colors [log(F 25 /im/T’i 2 /im) > 
0.5]. Unfortunately, however, the sensitivity of Charnley &: Latter (1997) was exceptionally 
shallow only for this object (the achieved rms was 0.0043 K in [K]). If OH 231.8-1-4.2 
exhibits the methanol line, of which the intensity is equivalent to that of 119312, the emission 
could be most likely missed due to the expected low signal-to-noise ratio of about 2. Therefore, 
we note that it is worth to make a more sensitive search toward OH 231.8-1-4.2 for conhrmation. 
We also note that the distance to OH 231.8-1-4.2 is somewhat closer to us(1.5 kpc; Choi et ah 
2012 ) compared to that to 119312 (3.8 kpc); if the 12 micron IRAS flux is normalized at 3.8 kpc, 
the flux of OH 231.8-1-4.2 is roughly 7.6 times weaker than that of 119312. 

In addition, some sensitive searches covering methanol lines have been undertaken toward 
the well-studied carbon star, IRC-l-10216 (Kawaguchi et ah 1995; Cernicharo et ah 2000; Ford 
et al. 2004; He et al. 2008), and the search results were negative in all cases. Even though 
the characteristics of IRC-l-10216 are somewhat different from 119312, we shall compare the 
sensitivities of these observations with the present case, because the methanol searches toward 
IRC-l-10216 are the most sensitive search toward evolved stars. 

Kawaguchi et al. (1995) made a line survey in the 30-50 GHz range using the Nobeyama 
45 m telescope; the observed frequency range included the CH 3 OH 1(0, 1)—0(0, 0) A^^ and 
7(0, 7)—6(1, 6 ) A^^ lines, which are detected in the present observation toward 119312. The 
rms noise achieved by Kawaguchi et al. (1995) was 0.008-0.011 K (T*). If we observed 119312 
at this rms level, the expected signal-to-noise ratio would be 4.1a - 5.6a for the 1(0, 1)—0(0, 
0) A++ line, and 17.7a - 24.4a for the 7(0, 7)—6(1, 6 ) A++ line. With this high signal-to-noise 
ratio expected, we can easily understand the peculiarity 119312 in its intensity of methanol 
emission. We should note that the distance to IRC-l-10216 is 120-150 pc (Groenewegen et al. 
1998; Lucas & Guelin 1999), while that to 119312 is 3.8 kpc (Imai et al. 2011): i.e., 119312 is 
even much further distant than IRC-l-10216. 

The frequency ranges of Cernicharo et al. (2000); Ford et al. (2004); He et al. (2008) 
(129 GHz-267.5 GHz) do not overlap with that of the present observations. However, since our 
model calculation predicts intensities in higher frequencies, we briefly compare the sensitivities 
about the high frequency lines. In table 7, we show a list of the methanol lines that are predicted 
to exhibit a strong intensity brighter than ~ 0.3 K; these lines are in the frequency range 
of Cernicharo et al. (2000); Ford et al. (2004) and He et al. (2008) except for the two lines at 


193 GHz. 

For the 145 GHz and 157 GHz lines [3(—1, 3)—2(—1, 2) E, 3(0, 3)—2(0, 2) A++, and 
1(0, 1)—1(—1, 1) E], the rms noise levels achieved by the previous observations at 145 GHz and 
157 GHz were respectively 0.0135 K and 0.0117 K (Gernicharo et ah 2000) and 0.008 K and 
0.031 K (He et ah 2008) both in Tmb- This means that if the model prediction is correct, the 
145 GHz and 157 GHz lines should be detected respectively at the signal-to-noise levels of 61a 

- 118(7 and 14(7 - 38(7. For the 241.8 GHz lines [i.e., 5(—1, 5)—4(—1, 4) E, and 5(0, 5)—4(0, 4) 
A"*"], the rms noise levels achieved by the previous observations were 0.002 K (Ford et ah 2004) 
and 0.015 K (He et ah 2008) both in T^b- Therefore, the expected signal-to-noise level is 22a 

- 198(7 for this line. 

As described above, we have to say that the properties of methanol lines in 119312 are 
far different from those expected in evolved stars. If 119312 is really an evolved star, the hidden 
mechanism, which causes the peculiarly strong methanol emission, is likely not directly related 
to the nature of the star, but rather the ambient materials and/or the interaction with ambient 
material presumably play a role. 

4-2. Characteristics of Methanol Emission 

The properties of detected methanol thermal lines are not consistent with those of a 
spherically expanding component, which is lying at the center of the nebulosity (Nakashima & 
Deguchi 2005; Murakawa et ah 2007, see, figure 3). There are two pieces of evidence supporting 
this statement. (1) The line prohles of the methanol thermal lines are not consistent with the 
kinematics of the spherically expanding component. Our previous radio observations of 119312 
(Nakashima et al. 2004; Nakashima (fe Deguchi 2005) revealed that the spherically expanding 
molecular component exhibits an expanding velocity of about 40 km s“^. This expanding 
velocity is clearly inconsistent with the narrow linewidth of the detected methanol thermal 
lines showing a linewidth of 3-4 km s“^. (2) The derived kinetic temperature of methanol 
lines (20-40 K) is most likely not consistent with the temperature of the spherically expanding 
component. This is because the dust temperature of the spherically expanding component is 
estimated to be 150 K-200 K using the mid-infrared photometric data (Nakashima et al. 2004). 
The high dust temperature presumably suggests that the spherically expanding component 
shows a much higher kinetic temperature than that measured using methanol lines. 

Based on the estimated kinetic temperature, we can roughly infer the location of the 
methanol emission region. The low temperature naturally suggests that the emission region 
should be at the outer region of the nebulosity, which is outer than the spherically expand¬ 
ing component. Our previous radio interferometric observation revealed that the spherically 
expanding component is surrounded by an outer expanding component with a size of about 
15"-20" (see, hgures 7 and 8 in Nakashima & Deguchi 2005; see, also hgure 3); the outer ex¬ 
panding component exhibits a bipolar structure. The expanding velocity of the outer expanding 
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component is about 3-5 km s“^, which is consistent with the linewidth of the methanol thermal 
lines. The kinetic temperatures of the outer expanding component, which are determined by 
the intensities of CO and NH 3 lines, are 13 K and 19 K, respectively (Deguchi et ah 2004); 
those are very close to the lowest possible kinetic temperature obtained from methanol lines. 
Therefore, the outer expanding component seems to be one of the possible sites of methanol 
molecules. 

The detection of class I methanol maser lines gives us a further hint about the site of 
methanol molecules. It is well known that there are two classes in methanol masers (i.e., class 
I and II), and the emission region of the two classes of masers (and therefore excitation mech¬ 
anisms) are different each other (see, e.g., Menten 1991). Class I masers are predominantly 
detected in shocked molecular gas (Kurtz et ah 2004), while class II masers are caused by radi¬ 
ation pumping. In the case of star forming regions, class I methanol masers are often detected 
along the periphery of a bipolar outflow, which induces shocked molecular gas (Voronkov et ah 
2006). According to these properties of methanol maser, there possibly exists a shocked region 
(or regions) in the system or 119312. 

We would like to note that even small relative velocity (less than 10 km s“^ ) between 
two colliding components can cause class I masers (Voronkov et ah 2006). In fact, in general, 
class I methanol maser emission is observed very close to the systemic velocity of the source. 
This suggests that in most cases the shocks are low velocity and/or primarily perpendicular to 
the line of sight. This fact suggests that there are three possible sites of class I methanol masers 
within the system of 119312. The hrst possible site of class I maser is, of course, the periphery 
of a bipolar flow (Nakashima et ah 2011). Even though the spherically expanding component 
does not seem to include methanol molecules as discussed above, the bipolar flow can generate 
masers if the tip of the jet goes beyond the outer edge of the spherically expanding component 
[see (a) in hgure 3]. Then, the periphery of the spherically expanding component also may be a 
possible emission region of class I methanol masers [see (b) in hgure 3], because the expanding 
velocity of the spherically expanding component is 20-30 km s“^ (Nakashima & Deguchi 2005), 
which seems to be sufficient enough to induce class I masers by a shock interaction. 

Another possible site of class I masers may be the periphery of the outer expanding 
component [see (c) in hgure 3], because we recently detected ambient molecular gas outside 
of the outer expanding component by an on-the-hy (OTF) mapping in the CO J = 1-0 line 
(Nakashima et al., in preparation). Even though the expanding velocity of the outer expanding 
component is relatively small (3-5 km s“^), shocked molecular gas still may be induced if the 
gas motion in the perpendicular direction is assumed. Interestingly, near-infrared images taken 
by the Subaru 8 m telescope and the University of Hawaii 2.2 m telescopes revealed that the 
infrared nebulosity shows a cometary tail extended in the south-west direction (Deguchi et al. 
2004; Murakawa et al. 2007); this may suggest that 119312 may have a relative velocity against 
to the ambient molecular gas along this direction. 
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The true origin of the outer expanding component of 119312 would depend on the loca¬ 
tion of the emission site of methanol masers. If the cases (a) and (b) are the case, methanol 
molecules must be preexisting in the outer expanding component. In such a case, the origin of 
the outer expanding component must different from that of the spherically expanding compo¬ 
nent and the bipolar jet, because it is hard to explain the existence of methanol molecules in 
a stellar envelope of a mass-losing evolved star. On the other hand, if class I maser is induced 
in the periphery of the outer expanding component (and methanol molecules are not included 
in the outer expanding component), the outer expanding component is still possibly explained 
with a common origin with the spherically expanding component and bipolar jet, at least, in 
terms of the methanol chemistry. Thus, a high angular-resolution radio interferometry, which 
enables us to specify the location of the emission regions of class I maser, would be helpful to 
constrain further the interpretations of 119312. 

Until now, class I methanol masers have been detected toward high-mass star formation 
regions (see, e.g. Kurtz et ah 2004), lower-mass star formation regions (see, e.g. Kalenskii et 
al. 2010), interaction between supernova remnant and molecular clouds (see, e.g. Pihlstrom 
et al. 2011) and starburst galaxies (see, e.g. Ellingsen et al. 2014). Even though there have 
been no conclusive evidences about the origin of 119312 so far, many observational properties 
of 119312 are consistent with those of an a mass-losing evolved star showing a bipolar flow (see, 
e.g., Nakashima et al. 2011). If the central star of 119312 is an evolved star, the class I maser 
toward 119312 should be classihed into a new category of class I methanol masers, which is 
caused by an interaction between the outflow of an evolved star and ambient materials. 

5. Summary 

We have reported the results of a systematic methanol observation of the enigmatic IRAS 
object 19312-1-1950. We detected 6 methanol lines including 4 thermal lines at 48.372 GHz, 
96.739 GHz, 96.741 GHz, and 96.745 GHz, and 2 maser lines at 44.069 GHz and 95.170 GHz. 
The detected maser lines belong to the category of class I masers, which are known as a tracer 
of shocked molecular gas. For better understanding of the physical conditions of the emission 
region, we constructed a simple radiative transfer model, and derived some physical parameters 
including kinetic temperature and relative abundance. The line prohles and derived excitation 
temperature suggest that the methanol emission region does not coincide with the spherically 
expanding molecular outflow lying at the center of the 119312 system. The detections of class 
I methanol masers suggest that there is a shocked gas in the system of 119312. We gave 
an implication about the location of possible methanol emission regions based on the present 
results. If the central star of 119312 is an evolved star, as suggested in the past, the class I 
maser detected in the present observation is the hrst case detected in an interaction between 
evolved star outflow and ambient molecular gas. 
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Appendix 1. Other Thermal Lines Detected in the Present Observation 

In the present observations, we detected 3 thermal lines other than methanol lines. 
In table A-1, the line parameters of the detected 3 thermal lines are summarized, and in 
hgure A- 1 , the spectra of the lines are given. Among the 3 lines, the SO N, J = 3, 2 - 2 , 1 line at 
109.252212 GHz and the HC 3 N 12-11 line at 109.173638 GHz were already detected by Deguchi 
et al. (2004), and they discussed chemistry of the molecular component by constructing a line 
radiative transfer model. The HG 3 N 4-3 F = 5-4 line at 36.392365 GHz is newly detected in 
the present observation. The linewidth of the SO N, J = 3, 2-2, 1 line (~9-10km s“^) is slightly 
wider than that of the HG 3 N lines (~6-7km s“^). 
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Table 1. Receivers Used in the Observation 


Freq. range 

(GHz) 

Receiver 

System temp. 

(K) 

20.0-25.0 

H22 (HEMT) 

120-180 

35.0-50.0 

S40 (SIS) 

180-230 

42.0-44.0 

H40 (HEMT) 

200-250 

80.0-116.0 

TIOOV/H (SIS) 

120-250 


Table 2. Spectrometers Used in the Observation 


Frequency channel 
Frequency coverage 
Vel. resolution at 22 GFlz 
Vel. resolution at 43 GFlz 
Vel. resolution at 43 GFlz 


AOS-W 

AOS-H 

2048 

2048 

250 MHz 

40 MHz 

3.4 km s"^ 

0.5 km s"^ 

1.7 km s"^ 

0.3 km s"^ 

0.7 km s“^ 

0.1 km s“^ 


Table 3. Parameters of Detected CH 3 OH Lines (temperatures are in T*). 


Rest freq. 

(GHz) 

Transision 

^^eak 

(km S“^) 

-^peak 

(K) 

dint 

(K km s“^) 

rms 

(K) 

Integ. time 

(sec) 

44.069476*1' 

7(0,7)-6(l,6) A++ 

36.1 

0.195 

0.560 

0.037 

560 

48.3724561 

l(0,l)-0(0,0) A++ 

35.8 

0.045 

0.249 

0.018 

2540 

95.1695161 

8(0,8)-7(l,7) A++ 

36.0 

0.021 

0.115 

0.007 

3040 

96.7393631 

2(-l,2)-l(-l,l)E 

35.8 

0.054 

0.205 

0.009 

4300 

96.7413771 

2(0,2)-l(0,l) A++ 

35.5 

0.076 

0.228 

0.009 

4300 

96.7445491 

2(0,2)-l(0,l)E 

39.1 

0.012 

0.064 

0.009 

4300 


* Separately observed on 2011 February 26. 

^Gonversion factor to Jy at these frequencies is 2.89 Jy/K. 
^Gonversion factor to Jy at these frequencies is 4.00 Jy/K. 
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Table 4. Noise levels and Integration Times of Non-Detections (temperatures are in T*) 


Rest freq. 

(GHz) 

Transition 

rms 

(K) 

Integ. time 

(sec) 

23.121024 

9(2,7)-10(l,10) A++ 

0.010 

3200 

23.444778 

10(l,9)-9(2,8) A— 

0.010 

3200 

36.169290 

4(-l,4)-3(0,3) E 

0.013 

1280 

44.069476 

7(0,7)-6(l,6) A++ 

0.013 

940 

48.247572 

l(0,l)-0(0,0) E i/t = l 

0.022 

2540 

48.376889 

l(0,l)-0(0,0) E 

0.015 

2540 

94.541806 

8(3,5)-9(2,7) E 

0.006 

3040 

95.169516 

8(0,8)-7(l,7) A++ 

0.006 

3040 

96.755507 

2(1,1)-1(1,0)E 

0.004 

4300 

108.893929 

0(0,0)-l(-l,l)E 

0.010 

4300 


Table 5. Detected lines taken into consideration in the LVG calculation. 


Frequency 

(GHz) 

Transition 

'T^mod 

(K) 

(rms) 

(K) 

r 

48.372456 

l(0,l)-0(0,0) A++ 

0.06 

0.06 (0.02) 

1.65 X 10-2 

96.739363 

2(-l,2)-l(-l,l)E 

0.13 

0.13 (0.01) 

6.91 X 10-2 

96.741377 

2(0,2)-l(0,l) A++ 

0.18 

0.18 (0.01) 

7.12 X 10-2 

96.744549 

2(0,2)-l(0,l)E 

0.03 

0.03 (0.01) 

4.70 X 10-2 


Table 6 . Lines used only for upper limit confirmation. 


Frequency 

(GHz) 

Transition 

'j^mod 

(K) 

(Upper limit) 

(K) 

T 

23.121024 

9(2,7)-10(l,10) A++ 

0.00 

0.06 

1.1 X 10-® 

23.444778 

10(l,9)-9(2,8) A— 

0.00 

0.05 

-8.7 X 10-® 

36.169290 

4(-l,4)-3(0,3) E 

0.04 

0.09 

-3.9 X 10-2 

44.069476 

7(0,7)-6(l,6) A++ 

0.00 

0.08 

-3.2 X 10-3 

48.247572 

l(0,l)-0(0,0) E i/t = l 

0.00 

0.06 

2.6 X 10-3 

48.376889 

l(0,l)-0(0,0) E 

0.01 

0.06 

2.0 X 10-2 

94.541806 

8(3,5)-9(2,7) E 

0.00 

0.06 

1.1 X 10-"^ 

95.169516 

8(0,8)-7(l,7) A++ 

0.00 

0.06 

-8.1 X 10-4 

96.755507 

2(1,1)-1(1,0)E 

0.01 

0.01 

1.1 X 10-2 

108.893929 

0(0,0)-l(-l,l)E 

0.07 

0.02 

9.1 X 10-2 
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Table 7. Bright Lines > 0.3 K) expected from the model. 


Frequency 

(GHz) 

Transition 

'Y'uiod 

^br 

(K) 

145.097443 

3(-l,3)-2(-l,2) E 

0.827 

145.103194 

3(0,3)-2(0,2) A++ 

0.945 

157.270818 

1(0,1)-1(-1,1)E 

0.441 

193.441610 

4(-l,4)-3(-l,3) E 

0.633 

193.454370 

4(0,4)-3(0,3) A++ 

0.673 

241.767247 

5(-l,5)-4(-l,4)E 

0.333 

241.791367 

5(0,5)-4(0,4) A+ 

0.354 


Table A—1. Parameters of Detected LinesEother than methanol (temperatures are in T*). 


Molecule 

Rest freq. 

(GHz) 

Transision 

l^eak 
(km S“^) 

-^peak 

(K) 

tnt 

(K km s“^) 

rms 

(K) 

Integ. time 

(sec) 

SO 

109.252212 

N,J = 3,2-2,1 

36.3 

0.060 

0.385 

0.007 

4300 

HC 3 N 

36.392365 

4-3 F = 5-4 

36.8 

0.080 

0.301 

0.014 

1280 

HC 3 N 

109.173638 

12-11 

36.5 

0.095 

0.274 

0.007 

4300 
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Fig. 1. Spectra of detected CH 3 OH lines. The velocity for the second top panel is calculated with respect 
to the 2(0, 2)-l(0, 1) A++ line (the middle of the three detected lines). The vertical lines in the top and 
bottom panels indicates the intensity peaks of the detected lines. 
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Fig. 2. Visual summary of the model calculation (see, Section 3 for details). The yellow plus-mark 
represents the best fit solutions. 
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'Spherically expanding 
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Fig. 3. Schematic view of the molecular components of IRAS 19312+1950. The panels (a), (b) and (c) 
indicate different possible emission regions of class I methanol masers; the red broken lines represent 
the possible emission regions, in which shocked gas is generated by hydrodynamical interaction between 
different kinematic components. The relative sizes between different components is arbitrary; the actual 
sizes of the spherically expanding component and the outer expanding component (in the long axis) are 
roughly 10" and 20", respectively (see, Nakashima & Deguchi 2005). The actual size of the bipolar outflow 
is unknown, because so far it has been observed only in maser lines (Nakashima et al. 2011). The position 
angle of the jet axis is known to be roughly —37°. 
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Fig. A—l. Spectra of detected thermal lines, which are other than methanol lines. 
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